clear all
close all
clc

y = inline('((1+eta)-1)./eta', 'eta');

eta = linspace(-10^-9,10^-9,1000)';

figure;
plot(eta,y(eta));

y_vera = 1;

err_ass = abs(y_vera - y(eta));
err_rel = err_ass./abs(y_vera);

figure;
plot(eta,err_rel);